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Abstract. The propagation trajectories of ultra-high-energy cosmic rays (UHECRs) are inevitably affected by Galactic 
magnetic field (GMF). Because of the inevitability, the importance of the studies of the propagation in GMF have increased 
to interpret the results of recent UHECR experiments. This article reviews the effects of GMF to the propagation and arrival 
directions of UHECRs and introduces recent studies to constrain UHECR sources. 

Keywords: Ultra-high-energy cosmic rays, Galactic magnetic field 
PACS: 98.62.Hr, 98.70.Sa 



INTRODUCTION 

The origin of ultra-high-energy cosmic rays (UHECRs) 
has been an open problem in astrophysics for about 50 
years, despite both theoretical and experimental efforts. 
A main difficulty to identify UHECR sources originates 
from the charge of UHECRs, which leads to the deflec- 
tion of their propagation trajectories in Galactic magnetic 
field (GMF) and intergalactic magnetic field (IGMF). 
These fields make UHECRs arriving at the Earth lose 
directional information on their sources contrary to neu- 
tral particles, e.g., photons and neutrinos. However, the 
deflection is expected to be small (several degree) at the 
highest energies lO^'' eV) even at an upper limit value 
of averaged IGME Bigmf/cIGMf'''^ < (lnG)(lMpc)'/2 
[1], if UHECRs are protons. Recent 7-ray observations 
imply ~ 10^'^ G in voids [2], although this estimation 
is controversial [3], and give conservative lower limits 
of IGMF in voids lO'^-'*^ G [4, 5]. Although lower val- 
ues of these constraints seems to be positive to identify 
UHECR sources directly, structured IGMF, i.e., magnetic 
field in clusters of galaxies and filamentary structure, can 
significantly contribute to the deflections, which includes 
large uncertainty in the modelling [6, 7, 8, 9, 10]. 

The propagation of UHECRs in GMF has also been 
studied from early days [11, 12]. In 1990s, the impor- 
tance of GMF was revisited [13, 14], which was stim- 
ulated by potential anisotropy in the arrival distribution 
of UHECRs [15]. In addition, since UHECRs detected 
at the Earth are inevitably affected by GMF, the effect 
of GMF is essential to consider what UHECR sources 
are. Moreover, after the identification of several UHECR 
sources, UHECRs can be used as a background source to 
constrain local GMF structures. 

The importance of considering the effect of GMF has 
increased recently. In order to interpret the correlation 
between UHECRs and extragalactic astrophysical ob- 



jects reported by the Pierre Auger Observatory (PAO) 
[ 1 6] and to constrain UHECR sources, the propagation of 
UHECRs in Galactic space should be taken into account 
because of its inevitability. The properties of the propa- 
gation highly depends on the composition of UHECRs, 
which is still controversial between the PAO [17] and 
High-Resolution Fly's Eye (HiRes) [18] [and Telescope 
AiTay (TA) (see Dr Tameda's talk in this workshop)], 
since the deflections of trajectories of UHECRs are sen- 
sitive to their charge (e.g., Ref. [14]). Also, in the case of 
transient sources, GMF produces significant time-delay 
of UHECRs compared to neutral particles emitted at the 
same time [19]. 

This article reviews the effects of GMF to the propa- 
gation and arrival directions of UHECRs and introduces 
recent studies to constrain UHECR sources. 



GMF MODELING 

Here, we give readers a minimal summary on GMF mod- 
elling. See Dr Han's talk in this workshop for recent 
progress on GMF observations and more details of mod- 
elling. 

GMF is modelled by two components; magnetic fields 
in the Galactic disk and in the Galactic halo. Note that 
the halo field usually means magnetic field far from the 
Galactic disk. The disk field has been modelled by a 
bisymmetric spiral structure (BS) [20], an axisymmetric 
spiral structure (AS) [13], and concentric rings with field 
reversals (e.g. Refs. [21, 22]). Which model can repro- 
duce observations globally better is still a controversial 
issue [22, 23], but the existence of (at least one) field 
reversals has been reported (e.g. [24]). There is not evi- 
dence that the directions of the field above and below the 
Galactic plane is different. 
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FIGURE 1. Examples of the projected trajectories of UHECRs (left: protons, right: irons) with lO'^'^ eV with the same arrival 
directions in different GMF models. This figure is originally from Ref.[25]. 
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In the halo field, exponential decay models of the spi- 
ral disk field in the direction perpendicular to the Galac- 
tic plane have been used [13, 26]. Faraday rotation mea- 
surements of extragalactic radio sources indicates the di- 
rections of magnetic fields parallel to the Galactic plane 
above and below the Galactic plane are opposite [27]. 
This anti-symmetric structure is theoretically supported 
by so-called AO mode in dynamo theories. So, some 
models treat GMF in the disk and GMF in the halo sep- 
arately [28]. Recent observations have confirmed this 
treatment and have revealed the scale height of electron 
density in the disk higher than in previous understanding 
[22, 23]. 

A dipole field, as a GMF component parallel to the 
Galactic disk, is also often taken into account [26, 29]. 
The dipole field is motivated by (i) a vertical field of 
0.2-0.3 IJ.G observed in the vicinity of the solar system 
[30], (ii) strong filamentary magnetic field perpendicular 
to the Galactic plane in the vicinity of the Galactic center 
[31] and (iii) a theoretical result of AO dynamo, but there 
is no clear evidence of the dipole field. In extragalactic 
edge-on galaxies, X-shaped structures of magnetic field 
above and below the disk have been observed, which 
implies galactic wind [32, 33]. The shape of the vertical 
component of GMF is still a controversial issue. 

In addition to these coherent components, the turbu- 
lent component of GMF, whose strength is 0.5-2.0 times 
as strong as the coherent components, is implied [34]. 
Only upper limits of the correlation length of the turbu- 
lent component (~ 100 pc) have been estimated due to 
a limited angular resolution of radio observations [34]. 
Although the turbulent component is not important for 
protons, it plays an significant role for heavy nuclei be- 
cause the Larmor radius of heavy nuclei in the jJ-G field 
approaches to the correlation length. 



PROPAGATION OF UHECRS IN GMF 

In order to investigate the propagation of UHECRs in 
GMF, a backtracking method has been often adopted. 
In this method, particles with the charge opposite to 
that of particles which we are interested in are injected 
from the Earth and their trajectories are calculated. Then, 
we can regard the trajectories of these fl«r/-particles as 
those of particles arriving from extragalactic space. In 
many cases, all the energy-loss processes of UHECRs 
are negligible because of the size of the Galaxy smaller 
than the energy-loss lengths of UHECRs. This method 
allows us to save much CPU time, since only the particles 
reaching the Earth are taken into account. This method 
can be applied to simulate the arrival distribution of 
UHECRs taking GMF and even sti-uctured IGMF into 
account [29, 8]. 

Since the solar system is embedded in the Galactic 
disk 8.5 kpc away from the Galactic center, magnetic 
field in the vicinity of the solar system mainly contributes 
to the total deflections of UHECRs in Galactic space. 
Fig. 1 shows the examples of the trajectories of UHE pro- 
tons (left) and irons (right) with the energy of 10^^'^ eV 
with the same arrival directions (£, b) = (0.0°, -44.8°). 
The GMF models adopted in this figure are BS and AS 
models for spiral field in the disk with two exponential 
decay scale in the direction perpendicular to the Galactic 
disk. The symbols S and A means symmetric and anti- 
symmetric field below and above the Galactic plane. The 
trajectories of protons seems almost to be straight lines 
by eye, but are deflected several degrees mainly in the 
Galactic disk. The field reversals of the disk field af- 
fect the total deflection of protons. In the BS models, the 
nearest field reversal is located at 0.5 kpc inside the solar 
system. When a charged particle passes through the field 
reversal, the direction of its deflection becomes opposite, 
and the total deflection angle is suppressed [26]. There- 



fore, BS models, in general, predict smaller deflection 
angles of UHECRs than AS models. We can observe this 
effect even for irons (green line), but sometimes heavy 
nuclei are confined in the Galactic disk even in the high- 
est energy range. The propagation trajectory of a parti- 
cle shown by the red line is very complex and is con- 
fined in the Galactic disk for a long time, and therefore 
directional information on its source is completely lost. 
The distribution of the deflection angles of UHE irons 
was well studied by Refs. [14, 37], similarly to proton 
cases [38], demonstrating that irons with the energy of 
6 X 10^^ eV are deflected more than 25° in more than 
80% region of the sky assuming a GMF model proposed 
by Ref. [28] plus weak dipole field. Ref. [25] calcu- 
lated the back-tracked positions of UHECRs observed by 
the Pierre Auger Observatory (PAO) onto the surface of 
Milky Way sphere on the assumption of pure iron com- 
position and demonstrated the correlation between UHE- 
CRs and nearby galaxies [16] is strongly disturbed. Fig. 2 
is the arrival directions of 27 UHECRs above 5.7 x 10^^ 
eV in the first public data of the PAO [35] {red), the back- 
tracked directions of them on the assumption of protons 
{green) and irons {blue). A BS-S model is adopted. The 
black points are galaxies in the IRAS catalog within 75 
Mpc [36]. While the back-tracked directions are close to 
the original arrival directions in the case of protons, the 
back-tracked directions of irons are far away from the 
original directions and no longer correlate with the dis- 
tribution of galaxies. Pure iron composition at the Earth 
is an extreme and unrealistic case even in the case when 
all the UHECRs emitted from sources are irons because 
photodisintegration with cosmic microwave background 
during propagation in intergalactic space. Nevertheless, 
it is an intriguing problem whether the correlation can be 
reproduced under heavy-nuclei dominated composition 
if the correlation is not a statistical fake. 

As mentioned above, the spiral field in the vicinity 
of the solar system mainly contribute to the total de- 
flections. Thus, reflecting the direction of local magnetic 
field, UHECRs from extragalactic sources in the north- 
ern Galactic hemisphere are deflected to the direction 
of Galactic south. This tendency can be seen in the fig- 
ure. In the case of anti-symmetry above and below the 
Galactic plane, the deflection directions are opposite. In 
both cases, since the Larmor radius of a charged parti- 
cle is proportional to its energy, the arrival directions of 
UHECRs from a source are arranged in the order of their 
energy if the coherent component of GMF is dominated 
[26, 29, 25]. 

The Liouville's theorem leads the isotropic distribu- 
tion of UHECRs at the Earth if UHECR flux is isotropic 
outside GMF. However, the distribution of the momen- 
tum directions of back-tracked UHECRs injected from 
the Earth isotropically is not isotropic because there are 
the trajectories of UHECRs which cannot reach the Earth 
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FIGURE 2. Arrival directions of the 27 PAO events [35] 
{red) and their back-tracked directions assuming protons 
{green) and irons {blue) in galactic coordinates. A BS model 
with symmetry above and below the Galactic plane is assumed. 
Black points are galaxies in the IRAS catalog within 75 Mpc 
[36]. This figure is originally from Ref. [25]. 



[26, 39, 38, 37]. In other words, there are directions in 
which sources can provide the Earth with UHECRs ef- 
ficiently (or inefficiently). Such an effect is called mag- 
netic lensing (de-lensing) [40, 41]. Interestingly, this ef- 
fect can occur even in turbulent magnetic field in a few 
% of the whole sky. 

SEARCHING FOR UHECR SOURCES 

The correlation between UHECRs and the positions of 
source candidates has been expected to give a hint to 
identify UHECR sources, and therefore a lot of efforts 
has been dedicated. Although the modifications of the 
arrival directions by GMF are not taken into account to 
avoid model dependence in many studies, several tests 
considered plausible GMF models and obtained positive 
results. In early days, Ref. [42] tested the correlation be- 
tween observed UHECR events and BL Lac objects tak- 
ing a BS model into account. The authors treated the 
charge of the UHECRs as a free parameter, and obtained 
a positive correlation signal in the case of +\ con- 
traty to no signal in the case of Q = — 1 , where Q is the 
charge of UHECRs in the unit of the absolute value of 
the electron charge. This implies that the composition of 
UHECRs is mainly protons. Recently, Ref. [43] exam- 
ined the correlation between the PAO events and y-ray 
sources detected by Fermi Large Area Telescope [44]. 
Although strong (^ 4(7) correlation was found, the mod- 
ification by GMF did not largely change the result be- 
cause a Cen A region is dominated in the signal. 

Ref. [45] tested the positional correlation between ob- 
served UHECR events and galaxies with extended ra- 
dio jets taking modifications of their arrival directions 
by several GMF models. The authors claimed that one 
third of UHECRs correlates with extended radio jets if 
they are light nuclei and proposed the idea that the re- 
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FIGURE 3. Histograms of the separation angles between 
200 simulated protons above 6 x lO'^ eV and their sources with 
10^"* Mpc^-' within 75 Mpc which gives maximal significance 
against random event distribution. The peak positions depend 
on the field reversals of GMF. This figure is originally from 
Ref.[25]. 



maining two thirds are heavy nuclei, i.e., a mixed com- 
position scenario. If radio galaxies are responsible for the 
observed UHECRs, a BS-S or BS-A model is supported. 
Interestingly, several trajectories of UHECRs pass close 
to Galactic magnetars and microquasars under the BS- 
S model. The propagation of UHECRs from Galactic 
sources was studied by Ref. [46]. UHECRs with lO'*^- 
10^^ eV from Galactic sources inside the solar system 
arrive from the northern Galactic sky and from the south- 
ern Galactic sky for a AS-S and BS-S models, respec- 
tively. In the case that the A-type parity of GMF is real- 
ized, the difference between the two models is not clear. 
These results can easily scale to the case of heavy nuclei 
at the highest energy by using the rigidity scaling. Thus, 
large-scale anisotropy is expected at the highest energy 
if heavy nuclei are dominated and the S-type parity of 
GMF is realized. 

Another approach is to estimate typical angular scale 
within which positive correlation between UHECRs and 
their sources appears by simulations. This angular scale 
can be used for the cross-check of the results from cor- 
relation studies between observed UHECRs and the dis- 
tribution of source candidates (galaxies), because we do 
not know what objects are UHECR sources. In addi- 
tion, since the typical angular scale depends on composi- 
tion, the estimation also gives an important information 
to constrain composition indirectly. Ref. [25] simulated 
the arrival distribution of protons taking GMF into ac- 
count and estimated typical angular scale within which 
positive correlation between the simulated protons and 
nearby sources used in the simulations. Fig. 3 shows one 
of the histograms of the angular scale within which the 
correlation signal is maximized in the case of the number 



density of UHECR sources of ris — 10^"* Mpc^^. Results 
a bit depend on GMF models and the number density 
of UHECR sources n^. The typical angular scale is ^ 5° 
for ris = 10"'* Mpc-3 and - 6° for n, = 10"^ Mpc^^ in 
the cases of BS models. AS models predict angular scale 
larger than those in BS models due to larger deflections 
in AS models. The smaller angular scale is predicted in 
the larger number density because the cross-correlation 
function includes correlation between simulated UHE- 
CRs and sources not emitting the simulated UHECRs lo- 
cated in the vicinity of the sources contributing to the 
observed (simulated) events. 

Motivated by ordered deflections by the regular com- 
ponent, the reconstruction of GMF structure along the 
line of sight and of the positions of sources has been stud- 
ied by simulations [47, 48]. A recent study concluded 
that about ten events from a source above 3 x 10^'* eV al- 
lowed to reconstruct the source position with an accuracy 
of 0.5° and the orthogonal component of the magnetic 
field with respect to the line of sight with an accuracy of 
0.6kpcZ-i [48]. 

Candidates of the arranged events were found in the 
data recently published by the PAO in a Cen A region 
[49] . The source position reconstructed by the arranged 
events is located at 8.5° from M87, a nearby powerful ra- 
dio galaxy. According to a statistical test, the probability 
that such a set of events are realized in random back- 
ground is '--^ 3 X 10^^. This scenario offers the possibility 
that the excess of UHECR events around Cen A is not 
produced by Cen A, and gave us not only indirect evi- 
dence that M87 is a nearby source of UHECRs but also 
a possible solution to the fact that highest energy events 
do not seem to correlate with the position of the Virgo 
cluster [50,51]. 

GMF produces the time-delay of UHECRs compared 
to neutral particles emitted at the same time. The time- 
delay plays an important role if the production of UHE- 
CRs is transient. For transient sources, this time-delay 
makes the apparent duration of UHECR bursts and we 
observe the sources as apparently steady sources in the 
observation time-scale of human beings [52, 19]. Com- 
bined by the estimation of the number density of UHECR 
sources on the assumption of steady sources [53, 54], the 
time-delay by GMF leads an upper limit of the rate of 
UHECR bursts as (60-3000) Gpc^^ yr"' [19]. Although 
the range means that this upper limit depends on the 
modelling of GMF, this limit rules out several transient 
phenomena for UHECR sources. 

SUMMARY 

We have reviewed the effects of GMF to the arrival direc- 
tions of UHECRs. The importance of GMF is in the fact 
that all the UHECRs arriving at the Earth are affected by 



GMF. A remarkable feature of GMF is its regular com- 
ponent. The regular component deflects the trajectories 
of UHECRs effectively and therefore significantly con- 
tributes to the arrival directions of UHECRs despite the 
smaller size of the Galaxy than the propagation distance 
of UHECRs in intergalactic space. The dominance of the 
regular component in UHECR deflections also allows us 
to use UHECRs as a spectrograph of GMF. Combined 
with increasing number of detected UHECR events and 
the improvement of the understanding of UHECR com- 
position, the studies of UHECR propagation in GMF will 
help us understand the origin of UHECR sources. 
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